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Abstract. We present new ATCA 17- and 24-GHz radio images and ESO-NTT optical spectra of the radio- 
loud Seyfert galaxy IC 5063, the first galaxy in which a fast (~ 600 kms" 1 ) outflow of neutral hydrogen was 
S_t ' discovered. The new radio data confirm the triple radio structure with a central, unresolved flat-spectrum core 

and two resolved radio lobes with steep spectral index. This implies that the previously detected fast outflow 
of neutral gas is occurring off-nucleus, near a radio lobe about 0.5 kpc from the core. The ionised gas shows 
highly complex kinematics in the region co-spatial with the radio emission. Broad and blueshifted (~ 500 kms" 1 ) 
emission is observed in the region of the radio lobe, at the same location as the blueshifted Hi absorption. The 
velocity of the ionised outflow is similar to the one found in H I. The first order correspondence between the radio 
' and optical properties suggests that the outflow is driven by the interaction between the radio jet and the ISM. 

However, despite the high outflow velocities, no evidence is found for the ionisation of the gas being due to fast 
t-H ■ shocks in the region of the outflow, indicating that photoionisation from the AGN is likely to be the dominant 

t-H ' ionisation mechanism. 

The outflow rate of the warm (ionised) gas is small compared to that of the cold gas, similar to what is found in 
I ■ other radio galaxies. The mass outflow rate associated with the Hi is in the same range as for "mild" starburst- 

' driven superwinds in ULIRGs. However, in IC 5063, the AGN-driven outflow appears to be limited to the inner 

kpc region of the galaxy. The kinetic power associated with the Hi outflow is a small fraction (a few xlO -4 ) of 
the Eddington luminosity of the galaxy but is a significant fraction (~ 0.1) of the nuclear bolometric luminosity. 
. | In IC 5063, the observed outflows may have sufficient kinetic power to have a significant impact on the evolution 

of the ISM in the host galaxy. 

Key words, galaxies: active - galaxies: individual: IC 5063 - galaxies: ISM 



1. Introduction 

Huge amounts of energy are produced through the accre- 
tion of material onto the super-massive black hole situ- 
ated in the centre of an Active Galactic Nucleus (AGN) . 
This energy is released into the surrounding medium in a 
number of different ways, ranging from collimated radio- 
plasma jets to UV and X-ray emission. The regions around 
an AGN are, therefore, highly complex and host a wealth 
of physical processes. Gas in different phases (atomic, 
molecular and ionised) is observed in this very hostile en- 
vironment. This gas, in particular through its kinematics 
and ionisation, can carry the signatures of the effect of the 
AGN on its surrounding medium. 
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* Based on observations with the ESO-NTT and with the 
Australia Telescope Compact Array 



The energy released from the nucleus can produce gas 
outflows with high velocities (thousands of kms -1 ). This 
has been observed in many AGN, from Seyfert galaxies to 
quasars. Gas outflows are detected as blueshifted absorp- 
tion or emission line wings in optical, UV and X-ray spec- 
tra, (see e.g. Crenshaw, Kraemer & George 2003, Krongold 
et al. 2003, Elvis et al. 2002, Holt et al. 2006 and references 
therein). More recently, fast and massive outflows of neu- 
tral hydrogen, detected as 21-cm Hi absorption against 
the central regions of radio-loud galaxies, have also been 
discovered (Morganti, Tadhunter & Oosterloo 2005a). 

All these different types of gaseous outflows are of 
great interest for a number of reasons. Understanding the 
driving mechanism(s) of the outflows is crucial for under- 
standing the physical mechanisms at work in the central 
regions of galaxies with an AGN. Furthermore, massive 
AGN-driven outflows are now suggested to dramatically 
affect the evolution of galaxies due to the large amounts 
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of energy they feed back into the interstellar medium (see 
e.g. Silk & Rees 1998, di Matteo et al. 2005). 

Gaseous outflows can be driven by super-winds asso- 
ciated with large starbursts, (Heckman 2002, Rupke et 
al. 2005a, b; Veilleux et al. 2005). However, in the case of 
galaxies with an AGN, radiation or wind pressure from 
the regions near the active super-massive black hole (i.e. 
a quasar wind) are the likely drivers of the gas outflows 
detected at X-ray and UV wavelengths. In radio-loud ob- 
jects, the interaction of the radio plasma with the (rich) 
gaseous medium in the direct vicinity of the active nu- 
cleus can provide another mechanism that can drive the 
outflows (see also discussion in e.g. Batcheldor et al. 2007). 

For the fast outflows of neutral hydrogen many un- 
resolved issues remain, in particular regarding their ori- 
gin and their quantitative effect on the ISM. To answer 
these questions, it is crucial to know where the outflows 
are occurring with respect to the AGN, and to know the 
characteristics of the gas in other phases (e.g. ionised). 
These H I outflows occur, in at least some cases, at kpc dis- 
tance from the nucleus (see Morganti et al. 2005a). This 
large distance from the AGN suggests that the interac- 
tion between the expanding radio jets and the gaseous 
medium enshrouding the central regions is the driving 
mechanism of such outflows. The associated mass outflow 
rates can be as large as ~ 50 Af© yr _1 , comparable (albeit 
at the lower end of the distribution) to the outflow rates 
found for starburst-driven superwinds in Ultra Luminous 
IR Galaxies (ULIRG). This would suggest that these jet- 
driven outflows of neutral gas in radio-loud AGN can in- 
deed have a significant impact - similar to the one of su- 
perwinds - on the evolution of the host galaxies. It is also 
important to be able to compare the characteristics of the 
H I outflows with those of the ionised gas. For example, 
in the radio galaxies 3C 305 and 3C 293 (Morganti et al. 
2005b, Emonts et al. 2005, Morganti et al. 2003) the mass 
of the gas in the ionised outflow is much less than that 
in the neutral outflow, making the neutral hydrogen the 
dominant and therefore the most influencing component 
in the feedback process. 




Fig. 2. The slit position of the NTT observations super- 
imposed on the HST image (from the public archive) ob- 
tained with WFPC2 through filter F606W (that includes 
bright emission lines like [Oin]). 



All these considerations motivated a new study of the 
radio-loud Seyfert galaxy IC 5063. IC 5063 (z = 0.0110) 
is classified as a Seyfert 2 and is classified as an early-type 
galaxy that shows a complicated system of dust-lanes (see 
Figure 1). In polarized light IC 5063 shows high polariza- 
tion in the near IR (Hough et al. 1987) and a strong, broad 
Ha emission (Inglis et al. 1993). Like in some other Seyfert 
2 galaxies this suggest that there is a broad-line region 
which is obscured from our direct view and the broad-line 
radiation is scattered into our line of sight by scatterers 
outside the obscuring regions. IC 5063 is among the most 
radio-loud Seyfert galaxies known, while its radio power is 
at the lower end of the distribution for radio galaxies. IC 
5063 is the first galaxy where a fast neutral outflow was 
discovered (Morganti et al. 1998). In the radio continuum, 
IC 5063 shows a triple structure of about 4 arcsec in size 
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Fig. 3. Vectors proportional to the 17 GHz fractional po- 
larisation (with a 3cr cut-off) superimposed on the contin- 
uum 17 GHz image. Vector position angle has been cor- 
rected for Faraday rotation. 

(about 1.3 kpcQ; see Morganti et al. 1998) aligned with 
the dust-lane. The neutral outflow occurs against the NW 
and strongest radio component (Oosterloo et al. 2000; see 
also below) . With new high-resolution and high-frequency 
radio observations obtained with the Australia Telescope 
Compact Array (ATCA) and deep long-slit optical spec- 
tra obtained with the ESO New- Technology Telescope 
(NTT), we aim to study the optical kinematics and ion- 
isation properties in relation to the radio structure. The 
kinematics (and the ionisation level) of the gas shed light 
on whether a strong jet/cloud interaction is responsible for 
the outflow or whether other processes must be invoked. 
This object also allows a detailed comparison between the 
properties of the ionised- and neutral-gas outflows. 

2. The radio structure: ATCA observations 

IC 5063 was observed with the ATCA using the 6-km con- 
figuration - the longest available for this radio telescope. 
The observations were carried out on 11 May and 8 July 
2004 at 17.8 and 24.8 GHz respectively. These high fre- 
quencies allow us to achieve a relatively high resolution. 
We took data simultaneously at 17.73 and 17.86 GHz in 
the first observation and at 24.77 and 24.90 GHz in the 
second, with a bandwidth of 128 MHz for each of these 
frequencies. These separate frequencies allowed us to im- 
prove (radially) the uv coverage. For both observations 
the total integration time on IC 5063 is about 9 h con- 
sisting of many cuts of 10 min spread over a 12-h period. 
PKS 1934-638 was used as the primary (flux) calibrator 
(with an assumed flux of 1.079 Jy at 17.8 GHz and 0.745 
Jy at 24.8 GHz). A nearby secondary calibrator (PKS 
2052-474) was observed for 1 minute every 10 minutes 
as well as for reference pointing. This reference pointing 
was done every hour. 



1 We assume Ho = 71 km s 1 Mpc 1 which implies a dis- 
tance of 46.5 Mpc; 1 arcsec is equivalent to 0.22 kpc 



The data reduction was done using the MIRIAD pack- 
age (Sault, Teuben & Write 1995). The final images were 
obtained using a few cycles of self-calibration and using 
uniform weighting to achieve the highest possible resolu- 
tion. The 17-GHz image has a beam size of 0.5 x 0.37" 
(12°) and an rms noise of 0.11 mJy beam -1 . At 24 GHz 
the beam is 0.36 x 0.29" (2.5°) with an rms noise of 0.18 
mJy beam -1 . 

Fig. 2 shows the 17- and 24-GHz continuum images. As 
expected from the previous 8-GHz observations (Morganti 
et al. 1998), the source shows a triple structure oriented 
along P.A.~ 115°. In contrast to previous lower resolution 
images, the bright NW source appears to be connected 
to the core by a low-brightness bridge. Although most 
of the flux in the NW source appears to come from a 
quite compact and bright region, low- brightness, diffuse 
emission surrounding the NW source is also detected. 

The new, higher-resolution observations allow us to in- 
vestigate in more detail the structure of the radio emis- 
sion while, due to the broad frequency range of the new 
observations, we are able to accurately derive the spectral 
index. The spectral index is crucial information for iden- 
tifying which of the features correspond to the core and 
which to radio lobes. 

The total fluxes of the three components at the two 
different frequencies are listed in Tablt[TJ The correspond- 
ing spectral indices (a defined as S = v~ a ) are also listed 
in this table. In contrast to the E and W component, the 
central component has a flat spectral index (a = 0.2), 
confirming the earlier prediction that this component is 
the radio core in IC 5063 (Oosterloo et al. 2000). The 
two east and west extended components show, instead, 
much steeper spectra, confirming their nature as radio 
lobes. Thus, the new observations support the interpreta- 
tion of the radio structure given by Morganti et al. (1998). 
Given that the neutral hydrogen outflow is detected (via 
VLBI observations, see Oosterloo et al. 2000) against the 
stronger (i.e. the western) radio component, the new ob- 
servations now confirm that this neutral outflow is occur- 
ring against a lobe located about 0.5 kpc from the core. 

We have derived the polarisation characteristics of the 
source at the two frequencies. We have obtained images of 
the Stokes parameters (Q, U), the polarised intensity (P = 
(Q 2 + C/ 2 ) 1 / 2 ) and the position-angle of the polarisation 
(x = 0.5arctan(£//Q)). The rms noise of the Q and U 
images is about 0.08 and 0.1 mJy beam -1 at 17 and 24 
GHz respectively. 

The polarised intensity and the fractional polarisation 
(m = P/I) were estimated only for those pixels for which 
P > 3<7qu- Polarised emission is found mainly over an 
elongated feature in the western radio lobe, stretching 
from the total intensity peak to the SW region of the lobe. 
Fractional polarisation values between 10 - 20% are seen 
towards the SW region at the two frequencies. The frac- 
tional polarisation decreases steadily towards the NE with 
the lowest values of under 1% measured in the centre of 
the lobe. Polarised flux is detected only in the western 
lobe. 
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Region 


17 GHz 


24 GHz 


a 




mjy 


mjy 




E 


5.6 


4.7 


0.53±0.02 


core 


30.8 


29.0 


0.18±0.02 


W 


91.2 


67.0 


0.93±0.02 



Table 1. Flux densities and spectral indices for the radio 
three components in IC 5063 



In Fig. 3 we have overlayed the position angle distribu- 
tion of the electric vectors corrected for Faraday rotation 
on the 17 GHz total intensity map. Using the values of the 
polarisation position angle at the two frequencies, we de- 
rive the rotation measure (RM) over the western lobe. We 
obtain high RM values over the region of up to a few thou- 
sand rad m~ 2 : at 5cr level the RM ranges between 2000 
and 3000 rad m -2 . We do not think the RM derivation 
suffers from nir ambiguities. At these high frequencies, the 
RM required for a rotation in position angle of tt radians 
is about 140000 rad m -2 - a very large and unlikely value. 

3. The ionised gas: NTT observations 

In order to obtain more information about the conditions 
of the gas in the inner regions of IC 5063, we have investi- 
gated the characteristics of the ionised gas using long-slit 
spectra obtained - in service mode - with the EMMI spec- 
trograph on the 3.5-m New Technology Telescope (NTT) 
in La Silla, Chile. Using EMMI in medium resolution 
mode, we obtained simultaneous spectra in the red and 
in the blue arm. The observations are summarised in Tab. 
2. In the blue, one central wavelength setting was used 
(A c = 4052 A) to obtain spectra covering 3560-4460 A. In 
the red arm, two settings for the central wavelength, A c 
= 5150 A and A c = 6315 A, were used yielding spectra 
with useful wavelength ranges of 4450-5760 A and 5700- 
7000 A respectively. Spectra were taken along PA 115° 
(see Figure [5]) with a 0.8-arcsec slit. This PA lies along 
the major axis of the radio emission as well as the dust- 
lane and includes both the nucleus and the radio hotspot, 
~2 arcsec to the west. The spatial pixel scale is 0.27 arc- 
sec/pix. To reduce the effects of differential refraction, all 
exposures were taken when IC 5063 was at low airmass 
(sec z < 0.1). 

The data were reduced in the usual way (bias sub- 
traction, flat fielding, cosmic ray removal, wavelength cal- 
ibration, flux calibration) using the standard packages in 
IRAF. The two-dimensional spectra were also corrected 
for spatial distortions of the CCD. The final wavelength 
calibration accuracy, calculated from the locations of the 
night sky emission lines (Osterbrock et al. 1996) is between 
0.05 and 0.1 A for the red spectra. The spectral resolution 
estimated from the widths of the night sky emission lines 
is 2.0 - 2.4 A for the red spectra. Due to the absence of sky 
lines in the blue spectra, we have not been able to quan- 
tify the wavelength calibration accuracy nor the spectral 



resolution in the blue, but there is no reason to believe 
that it is very different from that in the red spectra. 

Comparison of the data of several spectrophotometric 
standard stars taken with a wide slit (5 arcsec) throughout 
the run gives a relative flux calibration accuracy of about 5 
per cent. This accuracy is confirmed by the good matching 
of the continuum flux for the different regions. 

3.1. Line fitting and ionised gas kinematics 

Figure [4] shows the 2-dimcnsional spectra of IC 5063. Due 
to the very good spatial resolution and the high S/N, the 
2-dimensional spectra were analysed on a pixel-by-pixel, 
except in the outer regions where the lower S/N necessi- 
tated binning. 

Before proceeding, we would like to stress that in this 
paper we are mainly going to use the red spectrum and 
the analysis will focus on the kinematics of the emission 
lines. The main reason for this is that a proper analysis 
of the line ratios and diagnostic diagram would require 
a proper fitting of the underlying stellar continuum (as 
done in the case of both PKS 1345-12 and PKS 1549-79, 
Holt et al. 2003, 2006). This is important in order to take 
into account possible stellar absorption lines (if a young 
or intermediate stellar population is present, something 
that cannot be excluded for IC 5063 given the gas rich 
medium) . A first-order discussion of the emission lines ra- 
tios is given in Sec. 3.2 while a full analysis is postponed 
to a future paper. 

The spectra were extracted and analysed using the 
STARLINK packages FIGARO and DIPSO. Following 
the methods of Holt et al. (2003), we have modelled 
the emission lines, in each aperture, using multiple 
Gaussian components. We initially modelled the strong 
[O in] AA4959, 5007 doublet and defined a 'good fit' as that 
comprising the minimum number of Gaussian components 
to provide a physically viable good fit at each position. We 
will refer to these models as the '[Oiii] model' hereafter. 
The offset velocity is defined relative to the systemic ve- 
locity (~ 3400 kms" 1 ) derived from the Hi observations 
(Morganti et al. 1998). Up to four Gaussian components 
were required in the outflowing region while two or more 
components were required in other regions. Interestingly, 
the central region - corresponding to the AGN nucleus - is 
the least kinematically disturbed region. Examples of the 
[O in] emission line profiles in different regions are shown 
in Figure [5l The location of the radio lobe is also indi- 
cated for direct comparison between the characteristics of 
the ionised gas and the location of the radio emission. 

The spatial variations of the velocity widths and shifts 
and intensity of these components are shown in Figure [51 
We have attempted to group the various components us- 
ing their kinematical properties to aid in the description 
below. We realise that the separation in Gaussian com- 
ponents, as well as grouping them into "kinematical com- 
ponents" , can be somewhat artificial and subjective. We 
have based our grouping on kinematical continuity and/or 
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Fig. 4. The spectral regions from [On]AA3727 to [OlIl]A4363 (top), H/3A4860 to [N I]A5199 (middle) and 
[OI]AA6300,6363 to [S n]AA6716,6731 (bottom) along PA 115°. North-west is to the top, south-east to the bottom. 
Note, the spectra shown do not have identical scaling in either flux or spatial scale. Rather, the scales are optimised 
to highlight the differences in the intensity of the lines and allows the different kinematic components to be identified: 
from regular galaxy rotation at larger radii to the complex kinematics in the nuclear regions, clearly observed in all 
lines and is highlighted particularly well in the [O i]AA6300,6363 lines in this Figure. 
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Fig. 5. Spatial variations of the emission line profiles. In the centre is the 2D spectrum, zoomed in on the 
[O III] AA4959, 5007 (top) Ha/Nil blend (bottom). The four regions discussed in Section 4 are identified. On the left 
and right, a selection of the emission line profiles of [O in] AA4959, 5007 at various spatial positions are shown along 
with the best fitting model and the various sub-components required to model the lines. Key: black line (data), thick 
line (overall model) dashed lines are the different Gaussian used in the fit, see text for details. 
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Date 


Arm 


Exposure 


Set-up 


Slit PA 


Slit width 


A c 


AA 


Seeingf 


Conditions 


yyyy/mm/dd 




s 


CCD+grating 


degrees 


arcsec 


A 


A 


arcsec 




2000/07/25 


B 


4x1800 


TK1034+#12 


115 


0.8 


4052 


3560-4460 


0.8-1.0 


Photometric 


2000/07/25 


R 


2x1800 


TK2048+#7 


115 


0.8 


5150 


4450-5760 


0.8-1.0 


Photometric 


2000/07/25 


R 


2x1800 


TK2048+#7 


115 


0.8 


6351 


5700-7000 


0.8-1.0 


Photometric 



Table 2. Log of observations. Seeing is estimated using the DIMM seeing measurements and service observer's log. 
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Fig. 6. Spatial variations of the velocity width and shift of the various kinematic sub-components as measured from 
the [O in] AA4959, 5007 doublet. Crosses (red, in the electronic version) and open circles (green) represent the narrow 
components, filled triangles represent the extended broad component and open triangles (blue) represent the extremely 
broad component observed in the western hotspot. In addition, the points are scaled according to the measured line flux 
in the components. The solid horizontal line marks the position of the nucleus, and the dotted lines mark the position 
of the western (top) and eastern (bottom) hotspots. The systemic velocity used to estimate the "Offset velocity" is 
3400 kms" 1 as derived from Hi observations (Morganti et al. 1998). 



in similarities (e.g. FWHM) of the components. The kine- 
matics of the ionised gas are complex, in particular near 
the NW radio lobe. We identify the following main kine- 
matic components: 

1. As can be seen in Figures 2] and O the emission lines 
are highly extended along PA 115°. The brightest 
lines (Ha/ [Nil], [O m]AA4959,5007, [On]AA3727 and 
[Sn]AA6716/6731) are observed up to 15.1 arcsec (3.3 
kpc) to the E and 16.2 arcsec (3.6 kpc) to the W of 
the core. At these distances from the centre, a narrow 
component (FWHM ~ 100-200 kms" 1 ), which follows 
the rotation of the galaxy with velocity amplitude ~ 
200 kms" 1 is detected (the crosses in Fig. 6). 

2. In the region from the nucleus extending to at least 4 
arcsec SE and 2 arcsec NW, we detect a second nar- 
row component (FWHM ~ 200-300 kms" 1 , the open 
circles in Fig. 6) close to the systemic velocity. This 
component is seen outside the region of radio emission 
on the E side. 

3. A broad component, 500 < FWHM < 700 kms" 1 , 
is observed along the entire spatial extent of the ra- 
dio emission (the filled triangles in Fig. 6). In the 



region around the NW lobe (where we see the most 
extreme kinematics), this broad component appears 
highly blueshifted (~ 500 kms" 1 ), clearly deviating 
from the regular kinematics of the galaxy. 
4. A very broad component, (up to ~ 900-1300 kms" 1 ) 
only present around the region coincident with the NW 
radio lobe, and centred on the systemic velocity (i.e. 
not blueshifted, open triangles in Fig. 6) 

Hence, it is clear that the kinematics of the ionised 
gas is very complex. In particular, we note that there is a 
close correspondence between the ionised gas kinematics 
and the radio structure while the most extreme ionised gas 
kinematics occur, to first order, at the location where the 
fast outflow of neutral gas is detected. All these features 
suggest that a strong interaction between the radio plasma 
and the ISM is taking place. 

3.2. Physical conditions of the ionised gas 

A few interesting features can be extracted from the var- 
ious line ratios observed. We will concentrate on a few 
relevant regions (see Table [3]). 
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In the nuclear region, many strong emission lines 
are observed, as previously reported by Colina, Sparks 
& Macchetto (1991). We confirm the detection of high- 
ionisation species such as [Fex]A6375, [Fevn]A5721 and 
[Fe vn]A6087. In addition, we observe a feature at 5309A, 
similar to that observed in PKS 2152-699 by Tadhunter 
et al. (1989) which we identify as [Cav]A5309. The re- 
gion of the SE radio lobe is also characterised by strong, 
high-ionisation emission lines. In the region of the NW 
lobe, many bright emission lines are observed, although 
the highest-ionisation species found in the nucleus and 
eastern radio hotspot are not observed in this region. 

Optical images of IC 5063 show a complex dust lane 
system (see Figure^]) and, as our slit was aligned along this 
dust lane, it is likely that some of the emission line compo- 
nents will be reddened to some degree. We have estimated 
the degree of reddening in three regions (nucleus, north- 
ern outflowing region and the southern region with double 
narrow lines) using the commonly used Balmer decrement 
Ha/H/3. In all regions and in all kinematic components, we 
measure significant reddening, typically Ha/H/3 > 6, cor- 
responding to E(B — V) > 0.5. In all regions, the largest 
reddening is found in the component with intermediate 
blucshift and in one of the two narrow components (see 
Table [3]). Whilst our results are consistent with the pres- 
ence of a large, patchy dust lane, the specific numbers 
must be used with caution. No detailed modelling of the 
optical continuum was done and therefore no corrections 
for stellar absorption lines were applied to e.g. the H/3 
flux and it is therefore possible that we have underesti- 
mated the H/3 flux due to underlying stellar absorption, 
and so overestimated the degree of reddening. We have 
performed a quick fit to the stellar continuum to estimate 
the worst case scenario for possible underlying H/3 absorp- 
tion. In the southern and northern apertures, the effect of 
stellar absorption lines is negligible to a few percent and 
our estimated uncertainties on the emission line fluxes will 
cover this. In the northern region, the effect is again up 
to a few percent in the narrow components. However, the 
effect is perhaps as large as 20-25% in the intermediate 
component in the central region. Keeping this uncertainty 
in mind, we can make a few remarks about the line ratios. 

After modelling [O in] AA4959, 5007 in the various re- 
gions (see Sect. 3.1), we have attempted to model 
some of the other emission lines (limited to H/3 and 
[S n]6716/6731) using the same kinematic model as de- 
rived for the (Olll/ for the same aperture. We fixed the 
velocity widths and shifts and only allowed the relative 
fluxes of the kinematic subcomponents to vary. 

For one of the narrow components (the circles in Fig. 
6), the [O m] /H/3 ratio is high (~ 12) in the centre. As 
this component is narrow, it is unlikely to be strongly 
influenced by the interaction of the radio plasma with the 
ISM and the ionisation of this gas is likely due to the UV 
radiation from the nucleus. 

High [O m] /H/3 ratios (between 10 and 12) are also ob- 
served in the narrow components at the location of the ra- 
dio lobes. Furthermore, high values of the [Oiii]/H/3 ratio 



are also observed in the intermediate- width (600 kms -1 ), 
blueshifted component in the lobes. Even if no detailed 
analysis of the line ratios is done in this paper for IC 5063, 
the [O in] /H/3 ratios allow a useful comparison with the 
detailed studies of jet-cloud interactions in radio galax- 
ies. In the radio galaxies the [Oiii]/H/3 ratios provide the 
best evidence for shock ionization: the shocked compo- 
nents tend to have low [Oiii]/H/3 ratios (see e.g. 3C 171: 
Clark et al. 1997; PKS 2250-41: Villar-Martin et al. 1999), 
in contrast to the situation of the kinematically disturbed 
components to to the NW in IC 5063. Thus, this sug- 
gests that even if these components would originate from 
gas shocked by the interaction between the radio plasma 
and the ISM, these shocks are not the dominant ionisa- 
tion mechanism. Although this result will need to be con- 
firmed using more emission lines (and after a proper con- 
tinuum subtraction), it nevertheless supports the conclu- 
sions based on the energy budget (Morganti et al. 1998) 
and is also consistent with what has been found in similar 
studies of other Seyfert galaxies, e.g. Mrk 78 (Whittle et 
al. 2005). 

As with all other strong emission lines in IC 5063, 
the density diagnostic, the [S n]AA6716,6731 emission line 
doublet, is spatially resolved enabling us to investigate the 
density distribution of the emission-line gas. By modelling 
the doublet with the [O in] model, we derived the density 
sensitive ratio 6716/6731 A in the three regions mentioned 
earlier (Table 3). The narrow component seems to show 
a density of at most few hundred particles per cubic cen- 
timeter. For the NW region of the fast outflow, the best 
fit was obtained by setting the weak narrow components 
to the low density limit. The intermediate- width (~ 600 
kms -1 ), blueshifted component observed in the NW re- 
gion (as well as the intermediate component in the SE 
lobe) shows a much larger density (> 2000 cm -3 ). This 
is consistent with what is expected if these broad line are 
produced in the region of interaction between the radio 
plasma and the surrounding medium. 

4. Discussion 

The observations presented in this paper give strong ev- 
idence for the idea that the radio plasma ejected by the 
Seyfert nucleus in IC 5063 is strongly interacting with the 
ISM of this galaxy. The new, high-frequency radio obser- 
vations demonstrate that the brightest radio component 
is indeed a radio lobe. This implies that the fast, outflow 
of neutral gas is occurring off- nucleus, about 0.5 kpc from 
the core. The deep long-slit optical spectra show that the 
kinematics of the ionised gas is very complex over the en- 
tire region in which the radio emission is observed. The 
most extreme kinematics, in the form of a fast outflow of 
ionised gas, is detected at exactly the same location as 
the neutral outflow detected in Hi 21cm. On the other 
hand, despite the complex kinematics, the gas is most 
likely photoionised by the Seyfert nucleus. We now will 
discuss various aspects of these observations is some more 
detail. 
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Region component 



H/3 flux 

10 -16 



[Om]/H/3 Ha/H/3 [S ll]6716/6731 density E(B - V) 







erg s 1 cm 2 A 1 








cm 3 




SE 


narrow (x) 


7.8 ± 1.2 


9.3 ±0.2 


10.1 ±0.2 


1.04 ±0.05 


513±65 


0.94 




narrow (o) 


9.3 ± 1.2 


12.8 ±0.1 


11.1 ± 0.1 


1.06 ±0.05 


471±60 


1.02 




broad ( A ) 


11 (\ 4- 3 9 


y a it u. i 


0. 1 I U.o 


1.1 x u.uo 


qae;-i-i nn 
oyozt iuu 


U. / J 


centre 


narrow (o) 


31.0 ±2 


12.1 ±0.1 


12.9 ±0.1 


0.99 ±0.05 


632±80 


1.14 




narrow (x) 


2.8 ± 1.5 


4.6 ±0.5 


8.7 ±0.5 


0.99 ±0.19 


632+ 750 

uo ^-365 


0.82 




broad (a) 


(4.1 ±2.5) 


(31 ± 1) 


(26 ± 1) 


0.94 ±0.14 


774 +610 
' '^-340 


1.70 


NW 1 


narrow (x) 


0.9 ±0.6 


5.8 ±0.6 


1.0 ±0.6 


1.42 


< 100 


0.0 




narrow (o) 


1.7 ±0.4 


8.9 ±0.2 


4.5 ±0.3 


1.42 


< 100 


0.30 




broad (a) 


5.0 ±0.9 


9.7 ±0.2 


7.5 ± 0.2 


0.5 ±0.17 


> 3000 


0.70 




very broad (A) 


5.2 ±0.9 


6.9 ±0.2 


8.1 ± 0.2 


1.1 ±0.6 


395±52 


0.77 


NW 2 


narrow (x) 


2.8 ±0.2 


6.3 ±0.1 


3.5 ±0.1 


1.42 


< 100 


0.34 




narrow (o) 


1.2 ±0.2 


12.8 ±0.2 


7.6 ±0.2 


1.42 


< 100 


0.72 




broad (a) 


4.2 ±0.4 


13.3 ±0.1 


9.2 ±0.1 


0.7 ±0.1 


2390 ± 1000 


0.87 




very broad (A) 


6.4 ±0.4 


4.0 ±0.1 


6.1 ±0.1 


1.22 ±0.06 


213±25 


0.54 



Table 3. Emission lines measured in three main regions. Symbols are as in Fig. 6. The errors listed here are the 
combination of those associated with the fitting and due to flux calibration error. In bracket are very uncertain values. 
The fluxes for the SE region are derived from a 0.8 arcsec wide aperture. The NW represents two locations bracketing 
the hot spot. 



0.3 I r 




6650 6700 o 6750 6800 

Wavelength (A) 



Fig. 7. Example of fitting of the [S n] lines for the broad, 
blucshiftcd component 2 arcsec NW of the nucleus su- 
perimposed onto the data. The black line represents the 
overall model, the dotted line the narrow component, the 
dashed line the intermediate and the dot-dashed the broad 
component. 



4.1. Radio continuum and H I outflow 

Our new radio observations confirm the triple structure 
of the Seyfert galaxy IC 5063. From the estimated spec- 
tral index, it is clear that the central component corre- 
sponds to the radio core, while the other two compo- 
nents are radio lobes with a steep spectral index. The 
bright NW radio lobe is connected to the core by a faint 
bridge of radio emission that perhaps is a radio jet. The 
broad blucshifted H I absorption found in previous H I ob- 
servations (Morganti et al. 1998) is located - as deter- 



mined by VLBI observations - against this NW radio lobe 
(Oostcrloo et al. 2000). Our new radio data therefore con- 
firm that the outflow is occurring off-nucleus, i.e. against a 
radio lobe about 0.5 kpc from the nucleus. In this respect, 
IC 5063 is similar to other cases studied like 3C 305 and 
3C 293 (Morganti et al. 2005b, Morganti et al. 2003 and 
Emonts et al. 2005). 

The location of the fast outflow suggests that it is 
more likely caused by the interaction between the radio 
jet and the ISM than by other mechanisms such as e.g. 
radiation pressure. The interaction scenario assumes that 
a particularly rich medium is present around the NW 
radio lobe. That this is indeed the case is suggested by 
the detection, using NICMOS, of an asymmetric H 2 emis- 
sion, much stronger near the western lobe (Kulkarni et 
al. 1998) suggesting that an excess of molecular gas is 
present on the western side. Such emission is usually in- 
terpreted as evidence for fast shocks (for example, if the 
radio jet has struck a molecular cloud). The large rota- 
tion measures found in the W lobe also suggests that a 
particularly rich medium is present at the location of the 
western lobe. High values of the RM as found in IC 5063, 
are also measured in Compact Steep Spectrum sources. 
Twenty percent of the B3-VLA CSS sample of Fanti et 
al. (2004) have RM > 1000 rad/m 2 . In these cases, the 
Faraday screen most likely responsible for the RM is the 
magnetised interstellar medium. In most of these cases, 
the screen must be confined to the nuclear environment. 
An exception to this seems to be the high-redshift radio 
galaxy PKS B0529-549 (Broderick et al. 2006) where RM 
(-9600 rad/m 2 ) has been found. In this object, a gaseous 
halo is likely responsible for the for the extreme RM. The 
case of IC 5063 appears to be more similar to that of CSS 
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400 km/s 




sources. The presence of the high RM seems to further 
support the idea that the NW radio lobe is actually hit- 
ting a region of rich ISM. 

Using the column density of the neutral hydrogen 
(1 x 10 22 cm -2 assuming T spin = 1000 K) derived in 
Oostcrloo et al. (2000), we can estimate the mass of the 
H I involved in the outflow. Assuming that the absorption 
covers the NW radio lobe (~ 1 arcsec in size), we derive an 
H I mass of 3.6 x 10 6 M@. From the H I data the mass out- 
flow rate can also be derived (see Morganti et al. 2005a). 
Following Heckman (2002) and Rupke et al. (2002), the 
mass outflow rate for IC 5063 is estimated to be 35 Mq 
yr" 1 . As discussed in that paper, this is a substantial out- 
flow rate that is comparable (albeit at the lower end of 
the distribution) to that seen for starburst super-winds 
(Rupke et al. 2002) 



4.2. The ionised gas and the warm outflow 

The optical spectra of IC 5063 confirm the high ionisa- 
tion and the very complex kinematics of the ionised gas 
(Colina et al. 1991, Wagner et al. 1989). The region where 
the complex gas kinematics is observed is co-spatial with 
the radio emission. This is very similar to what was pre- 
viously found for other Seyfert galaxies and suggests an 
effect of the radio plasma on the surrounding medium, ei- 
ther via a direct interaction or via an expanding cocoon 
that surrounds the radio jet (e.g. Mrk 3, Capetti et al. 
1999). 

Interesting in the context of gas outflows, is the pres- 
ence of the highly blueshifted emission lines in the region 
of the bright NW lobe. We interpret this as a fast gas out- 
flow occurring in this region. Although we do not have a 
similarly detailed image of the complete region for the H I 
absorption, we know from the VLBI data that the most 
blueshifted Hi absorption component is also seen against 
the NW radio lobe. The similarity between the kinemat- 
ics of this blueshifted wing of ionised gas and the Hi is 
illustrated in Fig. [H Thus, a fast outflow is observed - at 
the same location and with similar velocities - in both the 
atomic neutral and the ionised gas. The fact that these 
outflows have such similar characteristics indicates that 
the detection of outflows of neutral hydrogen is likely not 
biased toward the strong radio continuum (i.e. one would 
see the outflow occurring only where there is strong back- 
ground continuum). Nevertheless, we can not entirely rule 
out the idea that the neutral outflow extends beyond the 
radio lobes and is driven by a different mechanism (e.g. a 
starburst-induced super- wind) . 

Following Holt et al. (2006), we can estimate the mass 
outflow rate of the warm, ionised gas and compare this 
with that of the neutral gas. We have determined the 
H/3 flux using a large aperture (1.5 arcsec) centred on 
the outflowing region to include all flux in the outflow- 
ing component. This gives a flux in the outflowing com- 
ponent of H/3 of (2.34 ± 0.27) x HT 15 erg s" 1 cm" 2 A" 1 . 
For the reddening and density, we have used values con- 
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Fig. 8. Comparison between the width of the Hi absorp- 
tion (white profile) and that of the ionised gas (from the 
[O ni]5007A). The first order similarity between the am- 
plitude of the blueshifted component is clearly seen. 



sistcnt with those presented in Table 3 for the outflowing 
region, namely a reddening of E(B — V) = 0.8 and a gas 
density of n e = 2000 cm~ 3 . This gives a reddening cor- 
rected H/3 flux of (3.53 ± 0.41) x 10" 14 erg s" 1 cm" 2 
A -1 and the luminosity of the outflowing component in 
H/3 (computed following the method of Holt et al. 2006) is 
L(H/3) = 9.77 x 10 39 erg s _1 cm~ 2 . A mass outflow rate of 
0.08 M Q yr _1 obtained. This is similar to the lower end of 
the range calculated for PKS 1549-79 (0.12-12 A/eyr" 1 ). 
Thus, regardless the uncertainties in the value of the den- 
sity and as earlier found in other galaxies, the contribution 
of the outflow of ionised gas is much lower than that of 
the Hi. As discussed in Morganti et al. (2005), the pres- 
ence of fast neutral outflows indicates that after a strong 
jet-cloud interaction the gas can cool very efficiently, as is 
indeed predicted by recent numerical simulations of jets 
impacting on gas clouds (Mellema et al. 2002, Fragile et 
al. 2004, Krause 2007). The relatively low outflow rate for 
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the ionised gas suggests that this cooling is indeed efficient 
and that most of the outflowing gas is neutral. 

It is intriguing to see that, while neutral and ionised 
gas share, to first order, the same location, only a rela- 
tively small part of the gas is ionised. Two possibilities 
can be considered to explain this: 

i) a large fraction of the jet-induced outflow is, for exam- 
ple, shielded from the AGN ionizing continuum. We could 
imagine a situation in which the jet drills through the rel- 
atively dense ISM we know to exist in the galaxy, staying 
relatively confined until it reaches a lower density part of 
the dust lane structure where the lobes/jets can expand 
more freely, causing the outflows we observe. Only a small 
fraction of the outflow is illuminated and ionized by the 
AGN continuum through the hole drilled out by the ini- 
tial expansion of the radio jet through the denser inner 
ISM. In this case, high resolution images of the ionized 
and neutral outflows would show that they have different 
detailed spatial distributions, while still both associated 
with the radio jets; 

ii) It is also possible that, following cooling behind the 
jet-induced shocks, the cooled cloud fragments are rela- 
tively large and dense. In this case, the emission lines could 
be emitted by the outer skins of the clouds photoioniscd 
by the AGN, with the masses/ volumes of the skins much 
smaller than the masses/volumes of the cloud fragments 
as a whole. 

Thus, in order to understand what is the situation we will 
need to image the exact location of the full Hi outflow, 
something that at the moment is not available. 

4.3. Impact of the outflow 

It is important to examine whether the outflows have a sig- 
nificant impact on the ISM of the host galaxy of IC 5063. 
Given that the H I outflow appears to be the dominant one, 
in terms of mass outflow rate, we will perform the calcula- 
tions using the parameters found for the neutral hydrogen. 
The kinetic power associated with the neutral outflow is 
of the order of 7 x 10 42 erg s _1 . This is derived by includ- 
ing the radial and the turbulent component of the flow, 
the latter estimated from the width of the line (see Holt 
et al. 2006 for more details). To investigate whether the 
characteristics of the observed outflows are interesting in 
the context of the feedback model (Fabian 1999, di Matteo 
et al. 2005), we can use the values derived from Nicastro, 
Martochia & Matt (2003) for the black-hole mass and ac- 
cretion rate of IC 5063. With a black-hole mass of 2.8 x 10 8 
M , the Eddington luminosity of IC 5063 is 3.8 x 10 46 erg 
s , this means that the kinetic power of the outflow rep- 
resents about few xl0~ 4 of the available accretion power. 
This result is similar to that found for PKS 1549-79 (Holt 
et al. 2006). However, unlike PKS 1549-79, which accretes 
at a high Eddington rate, Nicastro et al. (2003) found that 
IC 5063 accretes at a low rate (m ~ 0.02, defined as the 
ratio between the nuclear bolometric luminosity and the 
Eddington luminosity). Thus, in IC 5063 the kinetic power 



of the outflow appears to be a relative high fraction (<~ 8 
%) of the nuclear bolometric luminosity. Thus, the ratio 
between radiative luminosity (Lb i) and the kinetic power 
in the neutral wind appears to be consistent with that 
required in the Fabian (1999) model for the quasar feed- 
back effect. However, the overall power in the outflow is 
much less than envisaged in the numerical simulations of 
feedback in mergers, in which the quasar is accreting at 
close to the Eddington rate (e.g. di Matteo et al. 2005). 
Similar result is also found from the contribution of the 
warm absorbers as observed in the narrow line Syfcrt 1 
NGC 4051 (Krongol et al. 2007). 

4.4. The SEside 

What is the situation on the other side (i.e. on the SE 
side) of the nucleus? Also there an intermediate-width 
component, with velocity width ~ 600 km s- 1 is observed 
co-spatial with the radio emission. This suggests that, al- 
though not as extreme as on the NW side, also there the 
radio plasma interacts with the surrounding medium. As 
discussed above, the asymmetry seen in the conditions 
of the medium surrounding the radio source may explain 
some of the differences both in the radio morphology and 
in the kinematics of the g these asymmetries are 

likely to strongly affect the evolution of the radio plasma. 
We also see a component of ionised gas associated with 
the regularly rotating large-scale gas disk. In addition to 
all this, a second component is seen, blueshifted compared 
to the quiescent gas. It is hard to explain this component 
only by a jet-cloud interaction and it is perhaps more rem- 
iniscent of an expanding cocoon (see e.g. Mrk 3, Capetti 
et al. 1999). However, in the case of IC 5063, the second 
component is extending outside the radio emission and 
this makes the interpretation more complicated. We note 
that the large-scale gas disk in IC 5063 is seen edge-on, 
so the line-of-sight is intersecting the disk at various lo- 
cations, resulting in a large velocity widths and/or line 
splitting. This might explain the profiles at the SE side. 

5. Conclusions 

We have presented new radio and optical observations of 
the central regions of the Seyfert galaxy IC 5063. The ra- 
dio data confirm that the fast outflow of neutral hydrogen 
detected earlier, is indeed occurring off-nucleus, at about 
0.5 kpc from the core. The optical spectra show that the 
kinematics of the ionised gas deviates strongly from regu- 
lar rotation in the region coincident with the radio source. 
This supports the idea that the radio plasma ejected by 
the Seyfert core is interacting with the ISM of IC 5063. 
Very broad, blueshifted emissions lines are detected at the 
location of the NW radio lobe, at the same location at 
which the neutral outflow is detected, pinpointing the lo- 
cation of the strongest interaction. High rotation measures 
are detected in the NW lobe, suggesting the presence of a 
high density medium near the NW lobe. 
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Despite the strong kinematical evidence for an inter- 
action between the radio plasma and the ISM, we do not 
find, with the available limited data, evidence for the dom- 
inant ionisation being due to this interaction. Instead, 
photoionisation by the AGN is the most likely ionisation 
mechanism. 

While the mass outflow rate of the neutral gas is sub- 
stantial and is perhaps large enough for AGN feedback 
effects to be important in IC 5063, the mass outflow rate 
in the ionised gas is small compared to that seen in neutral 
hydrogen. This can be explained by very efficient cooling 
of the outflowing gas, as proposed in several models for 
jet-cloud interactions. 
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